Origin of the high-energy kink or the waterfall effect in the photoemission spectrum 
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The high-energy kink or the waterfall effect seen in the photoemission spectra of the cuprates is 
suggestive of the coupling of the quasiparticles to a high energy bosonic mode with implications for 
the mechanism of superconductivity. Recent experiments however indicate that this effect may be an 
artifact produced entirely by the matrix element effects, i.e. by the way the photoemitted electron 
couples to the incident photons in the emission process. In order to address this issue directly, we 
have carried out realistic computations of the photo-intensity in Bi2Sr2CaCu208 (Bi2212) where the 
effects of the matrix element are included together with those of the corrections to the self-energy 
resulting from electronic excitations. Our results demonstrate that while the photoemission matrix 
element plays an important role in shaping the spectra, the waterfall effect is a clear signature of 
the presence of strong coupling of quasiparticles to electronic excitations. 

PACS numbers: 79.60.-1, 74.25. Jb, 74.20.Mn, 74.72. Hs 



An anomalous 'high-energy kink' (HEK) in the disper- 
sion, which gives the associated angle-resolved photoe- 
mission spectrum (ARPES) the appearance of a 'water- 
fall' was first seen [1] in BizSrsCaCusOg (Bi2212) cuprate 
superconductor. Such HEKs or 'waterfall' effects have 
now been established as being a universal feature in the 
cuprates, and are found along both nodal and antinodal 
directions with their energy and momentum scales vary- 
ing with doping and details of the band structure [l|, Q- 
The HEK has been interpreted theoretically as providing 
evidence for interaction of the quasiparticles with some 
bosonic mode of the system0,0] with implications for the 
mechanism of high-temperature superconductivity. Re- 
cent experiments show however that the HEK is quite 
sensitive to matrix element (ME) effects, i.e. to the na- 
ture of the potoemission process or the way the incident 
photon couples with the electronic states of the system in 
generating the photoemitted electrons. In particular, the 
ARPES spectra undergo substantial changes in shape as 
one probes the electronic states by varying the energy of 
the incident photon or the momentum of the outgoing 
electron. In Bi2212, for example, the shape of the 

ARPES spectrum varies more or less periodically with 
photon energy from a shape that displays a single band- 
tail with relatively large intensity giving the spectrum a 
'Y' shape, to a spectrum which shows the presence of 
double-tails with a waterfall. [8| These results have led to 
speculation that the HEK may be an artifact produced 
entirely by ME effects.' 5, 8] On the other hand, for all 
polarizations and photon energies, the near Fermi sur- 
face features are strongly renormalized from LDA values, 
whereas the band bottom is not. Such a change in dis- 



persion with energy is generally taken to be a genuine 
signature of a bosonic mode coupling. 

In order to address this controversy, we have carried 
out extensive first principles, one-step photo-intensity 
computations in Bi2212 in which we include not only 
the effects of the ARPES matrix element, but also in- 
corporate a model self-energy based on accurate suscep- 
tibility calculations which properly reproduce the HEK 
phenomenology [3, 0, US]- In this way, we establish con- 
clusively that despite a strong modulation of the spectra 
due to the ARPES matrix element, a genuine HEK or a 
waterfall effect is still present in the cuprates, and that 
its presence indicates a significant coupling to bosons of 
electronic origin. Given the strength of the coupling and 
the high associated energy scale, these bosons are likely 
to play an important role in both the Mott and the super- 
conducting physics of the cuprates. In this connection, 
we also discuss model ARPES computations based on 
simplified tight-binding models for the purpose of gain- 
ing a handle on the interplay between the matrix element 
and self-energy effects, and also for delineating the na- 
ture of the striking characteristics of the ARPES matrix 
element such as the crossover from the Y-type spectral 
shape to the waterfall shape with photon energy. 

We discuss our key findings with reference to Fig. 1. 
Note first how the shape of the experimental spectrum 
changes dramatically at different photon energies. In 
panel (a) at 81 eV, the spectral intensity presents the 
appearance of a pair of 'waterfalls' with a region of low 
intensity through the middle of the figure. This is in 
sharp contrast to the measured spectrum in (b) at 64 eV 
where we see a 'Y-shape' with the two arms of the Y 
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FIG. 1: (Color online.) ARPES spectra in Bi2212. (a) and 
(b) are experimental results in the second Brillouin zone at 
photon energies of 81 eV and 64 eV, respectively.!^ (c) and 
(d) are corresponding theoretical photoemission spectra based 
on first principles, one-step calculations. 



connecting a vertical region of high intensity. Our realis- 
tic first-principles photo-intensity computations in which 
the matrix element as well as self-energy effects are ac- 
counted for reproduce the characteristic features of these 
shapes: The waterfall shape in panel (c) at 81 eV and 
the Y-shape in panel (d) at 64 eV, even though the arms 
of the Y in the computed spectrum in (d) are less clearly 
defined than the corresponding measurements in (b) . 

For the purpose of delineating the roles of the self- 
energy and matrix elements in shaping the ARPES spec- 
tra, Fig. 2 presents the results of photo-intensity com- 
putations based on a simplified two-band tight-binding 
Hamiltonian which models the low-energy electronic 
structure of Bi2212.[llilT^ The real and imaginary parts 
of the self-energy [l3| within the framework of the two- 
band model are shown in panel (a) for the antibonding 
band (AB) [similar results for the bonding band (BE) are 
not shown for brevity] . The model bands dressed by only 
the real part of the self-energy are shown in (b) , together 
with the LDA bands at — 2tt / c (magenta lines). The 
real part of the self-energy [solid blue line in (a)] is al- 
most linear in lu in the low-energy region. We write the 
slope of the linear part as (1 — Z~^) to define the renor- 
malization coefficient Z. This leads to a renormalized 
quasiparticle dispersion = ^^k, which, e.g., reduces 
the bi-layer splitting between the AB and BB bands at 
the antinodal point to Z{£^^q — £,^q), consistent with 
experimental results 0, [l3|. In contrast, at the F— point 
the dressed bands determined by ' 



move further away from each other. This is due to the 
fact that in Bi2212, ^^^/^^ hes in an energy range where 
the corresponding self-energy values S'(^p^^^^) possess 
an opposite sign. Moreover, the Green's function for the 
anti-bonding band shows multiple poles at the F— point 
because the line to — ^ coincides more than once with 
S'(a;) in (a), and as a result the anti-bonding band devel- 
ops an additional splitting. These opposing tendencies at 
low and high energies, which are seen clearly in the exper- 
imental spectra, are an unambiguous signature of strong 
coupling to a bosonic mode at intermediate energy. 
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FIG. 2: (Color online.) Model self-energy and spectral weight 
in Bi2212. (a) Real (blue solid line) and imaginary (red solid 
line) parts of the computed self-energy used in photo-intensity 
computations. The blue dashed line gives oj — ^ . (b) Dis- 
persion renormalized by the real part of the self-energy in 
(a) is compared to the bare dispersion of the anti-bonding 
band (AB, magenta solid line) and the bonding band (BB, 
magenta dashed line), (c) and (d) Spectral weights dressed 
by real and imaginary part of the self-energy for AB and BB, 
respectively, (e) and (f) Final photoemission intensities ob- 
tained after incorporating the matrix elements at two different 
photon energies, 95 eV and 75 eV, respectively. 

When the imaginary part of the self-energy S" 
is turned on, interesting spectral weight modulations 
emerge, which are shown separately in Fig. 2(c) and 2(d) 
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for the AB and BB bands, respectively. E" plays a cru- 
cial role in redistributing spectral weight such that the 
weight is shifted from the coherent region near the Fermi 
energy into incoherent parts at higher energies to pro- 
duce the HEK features seen experimentally [3|. However, 
the differences noted earlier near F lead to striking dif- 
ferences in the manifestation of the kink effect in the 
AB and BB spectra. In AB, the band bottom at F lies 
above the 'waterfall' region, and S" creates a long tail in 
the dispersion extending to high energies, so the fully 
dressed band exhibits the 'Y'-shaped pattern of panel 
(c). In contrast, the 'waterfall' shape emerges distinctly 
in the bonding band in panel (d) where the band bot- 
tom lies below the boson peak. These two dispersions 
are remarkably close to the two experimental dispersions 
observed at different photon energies as seen in panels 
1(a) and 1(b) Q and provide insight into the nature of 
the waterfall phenomenon in the spectra. 

In order to delineate contributions from different or- 
bitals, we have computed the photo-intensity within the 
framework of our tight-binding model using the linear 
electron-photon interaction. The matrix element for the 
Z^*'' band is then shown straightforwardly to be given 



X {i,lmn\fi)e' 



(1) 



Here kj is the momentum of the ejected electron, e de- 
notes the polarization of light with vector potential A, 
and Yim is the spherical harmonic for the angular vari- 
ables oikf. The final state is taken to be a free electron 
state. The initial state is a tight-binding state, which 
is expanded into atomic orbitals (nlm) of the i*'' atom in 
the unit cell at position Ri. The form factor 



Fniikf) - / r^drji{kfr)Rni{r), 



(2) 



where ji is a spherical Bessel function, is evalu- 
ated numerically using the radial part of the atomic 
wavefunction Il5l Il7l |. Fig. 3 shows that at low kj 
the oxygen contribution F21 for the oxygen and py or- 
bitals is dominant, while the Cu contribution F32 for the 
Cu dj.2_y2 orbital dominates for kf > 2 a.u.~^ . Although 
Eq. 2 is general, we have used only three orbitals, i.e. 
Cu d^2_y2, O Px^ and O Py in expanding the AB and BB 
bands to obtain the photo-intensities for our illustrative 
purposes. Finally, Eq. 1 can be recast into a useful form 
by collapsing all the symmetry information concerning 
the i*"^ orbital into the structure factor 5f (k/): 



M,(k^)=^5f(k;)e 



-ikf Ri 



(3) 
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3: Form factors of Eq. 1 for copper d^2_y2 orbital (solid 
F32, and for oxygen p orbitals, F2i. 



For a bi-layer system, the structure factor of Eq. 
becomes 



S^^S'c{h,k)[l± 



(4) 



where S'^ is the two-dimensional structure factor of Eq. 3, 
the -t- sign refers to the AB band and the - sign to the BB 
band. G = 27r{h/a, k/a, £/c) is a reciprocal lattice vector 
with a being the in-plane and c the out-of-plane lattice 
constant, d denotes the separation of the Cu02 layers in 
a bilayer. 18, 1^ The phase factor in Eq. 4 depends on 
the photon energy through 0, [20| 

kf = - Eund - $ + K)) - (fc|| + "||G||)2, (5) 

where hv is the incident photon energy and Ei,ind ~ 0.6 
eV is the binding energy of the electron in the solid at 
the waterfall, $ « 4 eV is the work function, Vb = 10 cV 
is the inner potential of the crystal and fcy -I- n\]^G\]^ is the 
total in-plane wave number, which we have taken to be 
« 27r/a to match the experimental conditions. 2l| Thus, 
whenever kfd changes by tt, the spectrum would switch 
from the odd to the even bilayer, a change that can be 
induced in view of Eq. 5 via the photon frequency v. 
This behavior is indeed seen in panels (e) and (f ) of Fig. 
2 where the matrix element is incorporated in the photo- 
intensity computations using Eqs. 1-5. In particular, at 
75 eV in panel (f), the AB band gets highlighted resulting 
in a Y-shaped spectrum with a tail extending to high 
energies. In contrast, in panel (e) at 95 eV, the bonding 
band dominates and spectral shape reverts to that of a 
waterfall with a double tail. 

Along the preceding lines. Fig. 4 further discusses the 
Y-waterfall spectral shape oscillation as a function of the 
photon energy. For this purpose, we consider in Fig. 4 
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FIG. 4: (Color online) Spectral weight integrated over the 
shaded binding energy window of Fig. 1(a) or 1(b) in the in- 
termediate energy region is shown to highlight how the spec- 
tra vary periodically between the Y- and waterfall shapes as a 
function of photon energy, (a) Experimental spectral weight 
normalized to the peak intensity at each energy. [§]; (b) The- 
oretical weights for the bilayer split spectra of Figs. 1(e) or 
(f). The colorbar is the same as Fig. 1. 

the integrated spectral weight over the shaded binding 
energy window shown in Fig. 1(a) or 1(b). The Y-shape 
is then characterized by a relatively narrow single tail in 
momentum (vertical axis in Fig. 4), while the waterfall 
displays a spHtting of this feature due to the presence of 
two tails. The experimental results of Ref. [sl shown in 
Fig. 4(a) are seen to be in good accord with the corre- 
sponding computations in Fig. 4(b), which are based on 
the bilayer dispersion of Ref. 22, some differences between 
theory and experiment with respect to the onset of Y- or 
waterfall shape in photon energy notwithstanding. No- 
tably, the theoretical results of Fig. 4(b) account for the 
effects of 'dimpling' of the Cu02 planes (i. e., deviations 
of the positions of the Cu and O atoms in the crystal 
structure of Bi2212 from a perfectly planar geometry) 
by including two bilayer splittings c?i and ^2 in model- 
ing the bilayer. When the self-energy correction is then 
included in the binding energy, Eund = Cfey + 5^(Cfc||), it 
makes fcj" complex in Eq. 5. This complex kj- damps the 
contribution from the lower Cu02 plane and thus modi- 
fies the intensity pattern through interference effects. We 
find that these corrections do not make much difference 
in our result. Taken as a whole, the comparisons be- 
tween the first-principles and model computations of the 
photo-intensity in Figs. 1, 2 and 4 also indicate that our 
simpler model with a limited number of orbitals is capa- 
ble of capturing many salient spectral features in the low 
energy region and giving insight into their origin. 

In summary, we have carried out extensive computa- 
tions of the photo-intensity in Bi2212 where the effects 
of the photoemission matrix element as well those of the 
coupling of the quasiparticles to electronic excitations are 
included realistically. We thus establish that despite the 
importance of the matrix element in shaping the spec- 



tra, the waterfall effect is a clear signature of the cou- 
pling of the electronic system to a high-energy bosonic 
mode, which bears on the physics of the pseudogap and 
the mechanism of high-temperature superconductivity. 
Our analysis of the spectra based on a simplified two- 
band tight-binding model reveals how the near-Fermi en- 
ergy bonding and anti-bonding bands associated with the 
Cu02 bilayers in Bi2212 produce characteristic Y-shape 
and waterfall shape of the spectrum as a function of the 
energy of the incident photons. Such a periodic modula- 
tion of the spectrum with photon energy may provide a 
new spectroscopic tool for getting a handle on the struc- 
tural aspects of the bilayer via the photoemission tech- 
nique. 
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